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ABSTRACT 

We study the physical and evolutionary properties of the “WD family” of AM CVn stars 
by computing realistic models of Interacting Double-Degenerate systems. We evaluate self- 
consistently both the mass transfer rate from the donor, as determined by gravitational wave 
emission and interaction with the binary companion, and the thermal response of the accretor 
to mass deposition. We find that, after the onset of mass transfer, all the considered systems 
undergo a strong non-dynamical He-flash. However, due to the compactness of these systems, 
the expanding accretors fill their Roche lobe very soon, thus preventing the efficient heating 
of the external layers of the accreted CO WDs. Moreover, due to the loss of matter from 
the systems, the orbital separations enlarge and mass transfer comes to a halt. The further 
evolution depends on the value of M after the donors fill again their lobe. On one hand, 
if the accretion rate, as determined by the actual value of (Mdon,Macc), is high enough, 
the accretors experience several He-flashes of decreasing strength and then quiescent He- 
burning sets in. Later on, since the mass transfer rate in IDD is a permanently decreasing 
function of time, accretors experience several recurrent strong flashes. On the other hand, for 
intermediate and low values of M the accretors enter directly the strong flashes accretion 
regime. As expected, in all the considered systems the last He-flash is the strongest one, 
even if the physical conditions suitable for a dynamical event are never attained. When the 
mass accretion rate decreases below (2 — 3) x 10“® Mq yr“^, the compressional heating of 
the He-shell becomes less efficient than the neutrino cooling, so that all the accretors in the 
considered systems evolve into massive degenerate objects. Our results suggest that SN .la or 
type la Supernovae due to Edge-Lit Detonation in the WD family of AM CVn stars should 
be much more rare than previously expected. 

Key words: Binaries: general, Supernovae:general, White Dwarfs, Accretion 


1 INTRODUCTION 

AM CVn stars are ultracompact cataclysmic binaries with spec¬ 
tra dominated by helium. At the time of writ ing 43 confirmed an d 
can didate objects were known, see Table 1 inlLevjtan et aP bOlSh . 
and IWagner et aljf20l4) : l^to. Hambsch & Monard bOl^ . Mea¬ 
sured orbital periods range from 5.5 to 65 min. There exist also 
several cataclysmic variables (some with hydrogen-deficient spec¬ 
tra) below the conventional minimum Porb o f CV (70-80 ) min. , 
which may be AM CVn stars in makin g llBreedt et all 201J; 
Carter_et_aflJ^ iLittlefield et ^ 1201.31 : iRamsav et alJ 2014 ; 
Gamavich et alj 2014h . The significance of AM CVn stars stems 
from their importance for the studies of very late stages of evo¬ 
lution of binary stars and of accretion disks physics; as well 
they are considered primary targets and verification sources for 


* E-mail: piersanti@oa-teramo.inaf.it (LP); tornambe@oa-teramo.inaf.it 
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high-frequency gravitational waves detectors. The current mod¬ 
els of AM CVn stars envision a semidetached binary harbouring 
a carbon-oxygen white dwarf (CO WD) accreting He-rich mat¬ 
ter. The donor may be either a helium WD or a low-mass helium 
star or a core of a main-sequence star strongly evolved prior to 
Roche-lobe overflow (RLOF). The evolution of AM CVn binaries 
is driven by angular momentum loss via gravitational waves radia¬ 
tion (GWR). An overview and a discussion of observational fea¬ 
tures, formation and evolution of thes e stars, as wel l as models 
for their disk s may be f ound, e. g.. iniw amej (|2 00^; iNelemans 
(2005, 2009|); Splheiml l20ip[); Ru itgLeLalJ l|2()l(j); Kottoet^ 


t20]2)\ Am^o^Seomieet^ r 2013 ); Postnov & YungelsonI ( l2014h : 
ICannizzo & NelemansI (l2015h . The topic of the present study are 
AM CVn systems with WD donors, sometimes also called “in¬ 
teracting double-degenerates” (IDD) or “white-dwarf family of 
AM CVn stars”. 

An essential issue defining the formation of AM CVn 
stars is stability of mass-transfer by degenerate donors 
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2 L. Piersanti et al. 


I Tutukov & Yungelsonl iNelemans et alj I 2 OOII : 

iMarsh. Nelemans & Steeghsll2004h . In double-degenerate systems, 
Roche lobe is first filled by the lighter of two WDs (in the context 
of this section - the secondary with mass M 2 ). Since for WDs 
the power of the mass-radius relation is close to —1/3, systems 
with initial q = M 2 /Mi > 2/3 are dynamically unstable. The 
components separation a of nascent AM CVn’s is so small that 
mass-exchange almost definitely begins in the direct-impact mode, 
without formation of the disk. Strict condition of stability for the 
case of no feedback of the angular momentum of the accreted 
material to the orbit becomes 


g < 1 -f (6 - ^i)/2 - + q)rc, 


( 1 ) 


where ^2 and /l are logarithmic derivatives of the radii of the donor 
and its Roche-lobe with respect to its mass, = Rc/a is the rel¬ 
ative circularisation radius. In between these two limits, dynamical 
stability of mass-transfer depends on the efficiency of spin-orbit 
coupling. 

An additional complication is brought in by the fact that the 
surface luminosity of the WD, as determined by the mass deposi¬ 
tion in the gravitational field of the accretor, can not exceed the 
Eddington limit. In the specific case of the gravitational poten¬ 
tial of interacting double-degenerates, if Ma is super-Eddington, 
the excess of the matter remains in the potential well of the 
accretor, the released energy heats it and may cause it to ex¬ 
pand and form a common envelope in which the two compo- 
ne nts will merge (I Han & Webbink|| l999l) . Numerical experiments 
by iMarsh. Nelemans & Steegh^ 2004 ) showed that in the case 
of weak tidal coupling stable mass transfer is possible if initial 
mass ratio is ^ 0.25 for Mco ~ 0.6 Mq and 0.21 for 
Mco « 1.0 M 0 . In the limit of very strong tidal coupling, the 
critical values of mass ratios become « 0.45 and « 0.3, respec¬ 
tively. 

The Eddington luminosity limit for M is more stringent than 
the dynamical stability one; therefore initial mass-transfer rates in 
AM CVn systems should be between 10“® and 10“® Mq yr“^. 
However, it is possible that, even if Ma < MEdd, it could ex¬ 
ceed the maximum rate of He-burning at the base of the He- 
envelope (lNomotolll982h . In the latter case a red-giant like ex¬ 
tended envelope forms and components, most probably, merg^B 
Thus, the existing interacting double-degenerates should start their 
evolution with M~ 10“® M© yr“^ and may evolve in a Hub- 
ble time to M^ 10~^^ M© yr~^ as it is theoretically inferred 
jTutukov & Yungelsonl 1 197!^ ; INelemans et alj 1200 ll ; IPelove et alj 
l2007h and sup ported by the analysis of the spectra of several 
AM CVn stars 1 Krausz et al.ll2010l ; fGehron et alj|2014l) . 


Along their path in Macc-M plane, accretors of IDD may ex¬ 
perience stable He-burning, burning via mild and strong flashes 
and, in principle, enter the regime of dynamical flashes (for a 
systematical study of the burning regim es of the accreted He see 
[Piersanti, Tornambe & Yungelsonll20l4 henceforth - Paper I). We 
define a flash “strong”, if in its course the WD oyerflows its Roche 


^ We note that a rigorous consideration of the net angular momentum 
chan ge following ballistic ejecti on of a particle, its motion and accre¬ 
tion iSepinskv & Kalogeral^l4h . as well as the accounting for angular 
mo mentum exchange between the spin s of the components and the or¬ 
bit (Kremer, Sepinskv & Kalogerdl2015h resulted in a slight change of the 
above-mentioned conditions for dynamical stability of mass transfer. In any 
case, the effects of the Eddington luminosity and the existence of a limiting 
He-burning rate are of higher value than these corrections. 


lobe. Otherwise, the flash is defined as “mild”. Interpolation be¬ 
tween the models computed with constant Ma suggests that in 
the course of the evolution AM CVn stars experience ~ 10 strong 
non- dynamical flashes. 

iBildsten et alj ( l2007tl paid attention to the circumstance that, 
with decreasing M, the mass to be accreted to get He-ignition, 
AMign, increases. Thus, there should exist the strongest “last 
flash”. Further, the mass of the donor becomes smaller than AMign. 
The last flash may become dynamical and result in a detonation, if 
the thermonuclear timescale, Tnuc = cpT/tnuc becomes shorter 
than the local dynamical time, Tdyn = Hp/cs, where Hp is the 
pressure scale-height and Ca is the sound speed. Due to the physical 
conditions existing in the accreting WD, the detonation produces 
short-living radioactive isotopes of Cr, Fe, and Ni. As well, ejected 
mass is small. The brightness of the event is comparable to sublu- 
minous SNe la, but its rise-time is only 2-10 day. Extrapolating the 
measured local birthrate of A M CVn stars and assu ming that all of 
them produce a visible event, IBildsten et alj ( l2007h estimated that 
the latter may occur once in 5000 - 15000 yr in a 10^^ M© E/SO- 
type galaxy., i. e., at a rate of ~ 1/lOth of the inferred SNe la rate 
in them. Having in mind the low lum inosity of thes e faint super¬ 
novae and the quoted occurrence rate. IBildsten et al] dubbed them 
“SN .la” (1/lOth of brightness at 1/lOth o f total rate). According to 
later calculations jShen & Bildstenl2009l) . such faint thermonuclear 
SNe are likely to occur if masses of accretors are in the range (0.8- 
1.2) M ©. We note that, accor ding to existing population synthesis 
results jNelemans et al.ll200ll) . AM CVn stars with such masses of 
the accretors should be extremely rare. Therefore, SNe .la may be 
substantially more rare events than estimated by Bildsten et al. At 
the moment, non e of several sugge sted SN .la candidates is defi¬ 
nitely confirmed dProut et al.ll2013l) . In the model of evolution of 
IDD above, it was assumed that mass-transfer is continuous. 

Recently, the authors of the present study suggested to account 
in the computation for the effect of mass and momentum loss from 
the binaries due to the Roche lobe overflow by the accretor dur¬ 
ing outbursts (Paper I). It was assumed that the matter leaving the 
system has specific angular momentum of accretofl 

In Paper I a grid of accretor masses and constant accretion 
rates was explored. In the present paper we consider the response 
to the accretion for time-dependent M, appropriate to AM CVn 
stars with WD-donors. As we show below, the RLOF episode and 
the associated mass and angular momentum loss from the system 
lead to the interruption of mass-transfer, thus resulting, after the 
flash, in epochs of cooling of accretor which change the character 
of thermal flashes and the course of evolution of the binaries under 
consideration. In Sect. we present our model and justify the se¬ 
lection of computed evolutionary sequences; results are described 
in Sect.[^ Discussion and conclusions follow in Sect.|4] 


2 SELECTION OF INITIAL BINARIES AND 
COMPUTATIONAL ASSUMPTIONS 

Considerations of the mass transfer stability, of the limits imposed 
onto the initial M by the Eddington luminosity and of the pos¬ 
sible formation of red-giant like envelope, suggest that the pre¬ 
cursors of WD-family of AM CVn stars are detached double- 


^ To some extent, a similar evolutionar y scheme was consi dered 
for semidetached W D+He-star systems by IVoon & LangeJ ^2004^ and 
iBrooks et al.l l2015l) : as well, it resembles the scenario of evolution of cata¬ 
clysmic variables with “hibernation”. 
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degenerate (WD+WD) systems with rather extreme mass ratios and 
massive accretors. Such binaries are observed as detached binary 
“extremely low-mass” white dwarfs (ELM), which harbour a He 
WD with mass below ~ 0.2 Mq and a much more massive (pre¬ 
sumably) CO companion; some of them have merger time less t han 
the Hubble time (see, e. g., iBrown et alj|2013l; iKilic et alJl2014l) . 

Currently, several systems are observed, which may be consid¬ 
ered as proper candidate precursors of AM CVn stars (note, here 
we call “primary” the more massive invisible component with Mi). 
The first candidate is eclipsing binary SDSS J075141.18-014120.9 
withPor-i,= 1.9hr, (M 2 -bMi)/M© = ( '0.19±0.02)-b(0 .97+° m ) 
and expected merger time ~160 Myr jKilic et ^ [20ll . Another 
“best” candidate is the eclipsing detached binary NLTT 11748 with 
Porb= 5.6hr and estimated mass of the secondary ranging from 
M 2 = (0.136 ± 0.007) to (0.162 ± 0.007) M© and, correspond¬ 
ingly, mass of the primary from Mi = (0.707±0.008) to (0.740± 
0.008)M© dKaplan et al.l l2014h . Several candidate systems are 
single-lined and for them only lower limits of Mi are available. 
Among them are SDSS J1741-I-6526 - Porb- 1.47 hr, M 2 = 
0.17Mo, Ml ^ 1.11 M© jKilic et al.ll2^ : SDSS J1141-1-3850 
- Port « 6.23 hr, M 2 = 0.17 M©, Mi ^ 0.76 M©; 
SDSS J1238-H1946 - Port ~ 5.46 hr, M 2 = 0.17Mo, Mi > 
0.64Mo; SDSS J 2132-t0754 - Port ~ 6.01 hr, M 2 = O.17M0, 
Ml > 0.95 M© jBrown et~^l2013h . Yet another candidate with 
less certainly estimated parameters is S DSS 1257-t5428 - Porb= 
4.55 hr, M 2 0.2 M© , Mi IM© jKulkami & van KerkwiitJ 


nr, IVI 2 u.zivi© , 
|2010|;lMarshetal.l201lL 


. At the end, other candidate AM CVn sys¬ 
tems may be hidden among observed sdB stars, deemed to evolve 
into CO WDs with mass 0.5 M©, with low-mass WD compan¬ 
ions, like PG 1043-1-760 - Port = 2.88 hr, M 2 > 0.101,M©, 
SPSS J083006-H47 510 - Port = 3.552 hr, M 2 ^ 0.137,M© 
jKupfer et alJ[2oT^ . 

Having in mind the parameters of candidate AM CVn sys¬ 
tems and for the sake of comparison with computations in Paper I, 
in the current work we selected the following donor and accretor 
combinations: (Mdon, Mace)=(0.17,0.60), (0.15,0.92), (0.20,1.02) 
M©. In the following, we address the sequences of models for these 
systems as S060-I-017, S092-I-015, and S102-I-020, respectively. The 
initial CO WD models have been obtained by evolving the “heated 
Models” M60, M092 and M102 from Paper I along the cooling se¬ 
quence down to the luminosity level L=0.01 L©. 

Some relevant physical quantities for these binary WDs are 
listed in Table [T] The donors are modele d as zero-temperature 
objec ts with Eggleton’s M-R relation (see IVerbunt & Rappaportl 
Il988[) . In the computation we neglect possible effects of finite ini¬ 
tial en tropy of the He WDs jPelove & Bildstenl2003l : lDelove et al] 
l2007h . but we discuss them in §4.. We do not consider the effects re¬ 
lated to the possible presence of a thin hydrogen la yer at the surface 


20061: Delove etalJ 2007 

1; iKanlan. Bildsten & Steinfadtj 

Shen. Guillochon & FolevI 

2 OI 3 I; IShenll2015l) and consider 


tion of He only. 


3 RESULTS 

At the beginning of the computation, the two components of 
each system are put in contact (i. e. Rdon*'*’ = Rdon). At each 
time step, the angular momentum loss via GWR is computed, so 
that the donor overfills its own Roche lobe. As a consequence, 
mass is removed in order to restore the condition Rdon = Rdon*'® 
and it is transferred conservatively to the accretor. If and when 


Table 1. Physical properties of the initial binary systems considered in the 
present work We list the separation a, the mass of the donor Mdon. the mass 
of the accretor Mace and some physical properties of the accretors, namely 
the mass extension of the He-deprived core Mco. the mass extension of the 
more external layer where the helium abundance by mass fraction is larger 
than 0.05 AMhc, the temperature Tc in K and the density pc in g cm~® 
at the center, the surface luminosity L, the effective temperature Tefj in K 
and the surface radius R. In the last row we report the cooling time tcool 
defined as the time elapsed from the bluest point along the loop in the HR 
diagram and the epoch of the initial model considered in the current work. 


LABEL 

S060+017 

S092-t015 

SI02+020 

a (in 10“^R©) 

7.599 

9.222 

7.965 

Mdon (in Mq) 

0.170 

0.150 

0.200 

Mace Mq) 

0.59678 

0.91962 

1.020467 

Mco (in M©) 

0.5158 

0.9114 

1.0158 

AMne (in IQ-^M©) 

1.551 

0.389 

0.056 

log(Tc) 

7.2820 

7.1478 

7.2329 

log(Pc) 

6.5626 

7.3378 

7.6109 

log(L/L©) 

-2.0077 

-2.1163 

-2.0322 

log(Teff) 

4.2102 

4.2640 

4.3148 

log(R/R©) 

-1.9012 

-2.0633 

-2.1228 

tcool (inl0®yr) 

1.647 

2.767 

2.889 


the accretor overfills its own Roche lobe, due to evolutionary 
reasons (e. g. the onset of a He-flash), mass is removed from 
the accretor, assuming that it is ejected from the system. The 
specific angular momentum of the lost matter is assumed to be 
equal to the orbital one of the accretor. The thermal response 
of the accretor to mass deposition is computed in detail by 
using an updated versio n of the FRANEC code, the original 
one being described in IChieffi & Stranier3 jl989ll . The setup 
of the code as well as the input physics are the same as in 
Paper I. The chemical composition of the matter transferred from 
the donor to the accretor is fixed as in Paper I, by assuming 
that all CNO elements have been converted into ^'^N, namely: 

Yi2© + Y'i 3 q -|- Y'i 4 j,j -f YI5J.4 ”'*^180 T Yisq = Y) 4 j^, 

where Yj is the abundance by number of the y'-isotope and the 
superscripts i and / refer to the initial MS star and the final 
He-donor WD, respectively. In Fig. [T] we show the profiles in the 
p — T plane for the accretors in our models at the beginning of 
the mass transfer process while in Fig. we plot as a function of 
the WD mass fraction the mass fraction abundances of "^He, 
and in the most external layers of the accretors for the same 
structures. 

When the mass transfer starts, the evolution of the accretors in 
the three considered systems is quite similar, as illustrated in Fig.[3 
where we show the evolutionary tracks in the HR diagram for the 
CO WDs in the investigated systems. Along the tracks, we mark 
with different symbols several important epochs. For each point we 
list in Tablel^some relevant physical properties of the accretors and 
of the binaries. 

When the matter falls onto the accretor, it delivers gravita¬ 
tional energy that is locally stored as thermal energy and triggers 
the evolution backward along the cooling sequence. The heating 
determined by mass deposition drives to He-ignitiorfl The max¬ 
imum luminosity attained during this phase depends mainly on 
the relative efficiency of the heating by accretion and the inward 


® Consistently with Paper I, we individuate the ignition epoch as the one 
when the nuclear energy per unit of time delivered by He-burning is 100 
times the surface luminosity. 
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0.1 1 10 100 1000 
p(10'‘ g cm“^) 


Figure 1. Profiles in the p — T plane for the accreting CO WDs at the be¬ 
ginning of the mass transfer process. Each panel refers to a different initial 
binary system, as labelled. Filled circles mark the He/CO interface while 
open ones the point where X(He)=0.05, roughly corresponding to the He- 
buming shell. 


M060 M092 M102 



M/M„o M/Mwd M/Mwd 


Figure 2. Mass fraction abundances of (solid lines), (dashed 
lines), and (dotted lines) in the most external layers of the accretors 
in the three considered binary systems. Abscissa represents the mass frac¬ 
tion. 



Figure 3. Evolution in the HR diagram of the accretors in the three binary 
systems listed in Table [J Along the tracks some points corresponding to 
certain epochs in evolution are marked by different symbols and letters, see 
Table[2|“FDC” means flash-driven convection. 


thermal diffusion; the former is determined hy the accretion rate, 
while the latter depends on the thermal structure of the He-buffer. 
As the He-burning shell has partially degenerate physical condi¬ 
tions, a thermonuclear runaway occurs, driving to a very powerful 
non-dynamical He-flash. As discussed in Paper I, this first flash 
represents a sort of “heating” mechanism which alters the thermal 
content of the whole He-buffer above the CO core and sets in the 
physical conditions suitable for quiescent He-burning. The latter 
could occur only after the flashing structures have attained their lo¬ 
cations in the HR-diagram corresponding to a post AGB-star with 
the same CO core and He-buffer masses. So this first He-flash rep¬ 
resents the analog in real binary systems of the “heating procedure” 
adopted in Paper 1. 

The energy delivered by the He-flash can not be removed via 
radiative transfer and, hence, a convective shell forms very soon 
(point A in the HR diagram) and rapidly attains the stellar surface 
(point D). Hence, the huge amount of nuclear energy delivered by 
the He-flash is injected into the whole He-rich buffer above the 
He-burning shell, whose thermal content becomes too large for a 
compact configuration as that of the flashing objects. In order to 
dissipate the too large thermal energy stored in the He-rich mantle 
the accretors start to expand. If the flashing structure would be iso¬ 
lated in the space, it could evolve freely, increasing its luminosity 
and, then, should expand to very large dimensions, thus dissipat¬ 
ing a part of its thermal energy via mechanical work and finally 
reaching at the He-shell the physical conditions suitable for qui¬ 
escent He-burning. But the objects considered in the present work 
are embedded in compact binary systems, so that they overfill their 
own Roche lobe very soon (point E). The following evolution of 
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Table 2. Physical properties of the accretors in the binary systems listed in Table^at same selected epochs during the first He-flash episode. The various times 
are: A - onset of the flash-driven convective shell; B - ignition of He-burning; C - maximum luminosity of the He-burning shell; D - flash-driven convective zone 
attains the stellar surface; E - beginning of the RLOF episode; F - flash-driven convective shell recedes from the stellar surface; G - flash-driven convective shell 
disappears; H - end of the RLOF episode; I - resumption of the mass transfer from the donor. For each epoch of each binary system we list surface luminosity 
L/Lq and temperature Teff, the mass coordinate of the He-buming shell Mhb iii Mq, the temperature Tne Ih® density pHe ■ cm~^ and 

luminosity Lhc in erg • of the He-buming shell, the separation of the system a in 10“^R©, the total mass of the accretor Mace in M©, and of the donor 
Mdon in M©, the time elapsed between two successive epochs At in yr. For each considered binary system we report also the amount of mass transferred 
from the donor to the accretor AMtran in 10~^M©, the mass lost during the RLOF AMiost in 10“^ M© and the corresponding value of the retention 
efficiency pace - Note, negative At for events D and F in S102-1-020 sequence mean that they happen before event E. The last block of data refers to the “Cold” 
S060-I-017 model described in Sect. l3.4l 



A 

B 

c 

D 

E 

F 

G 

H 

I 





S060+017 






log(Teft) 

5.128 

5.127 

4.873 

4.625 

4.560 

4.843 

5.311 

5.439 

4.864 

log(L/L 0 ) 

1.896 

1.893 

1.199 

0.454 

0.478 

1.622 

3.513 

4.024 

0.715 

Mne/M© 

0.60176 

0.60074 

0.59945 

0.59944 

0.59942 

0.59920 

0.59801 

0.59748 

0.58997 

PHe 

7.878 

7.896 

2.128 

1.250 

0.770 

0.174 

0.172 

0.255 

5.249 

Thb 

1.235 

1.343 

3.260 

3.539 

3.579 

2.893 

2.183 

1.638 

0.605 

log(LHe/L©) 

2.984 

3.894 

10.278 

9.745 

9.210 

6.406 

4.794 

2.738 

-6.119 

a (in lO^^R©) 

8.592 

8.592 

8.593 

8.593 

8.593 

8.734 

8.8230 

8.886 

8.504 

Macc/M© 

0.62559 

0.62560 

0.62561 

0.62561 

0.62560 

0.61316 

0.60557 

0.60018 

0.60018 

Mdon/M© 

0.14119 

0.14118 

0.14117 

0.14117 

0.14117 

0.14117 

0.14117 

0.14117 

0.14117 

At (yr) 

55637.0 

42.4986 

12.9458 

0.0006 

0.0018 

0.3910 

15.2389 

54.1418 

23945.8 


AMtran = 28.824 

AMio,t = 25.427 

Pace- 

=0.118 






S092+015 






log(Teft) 

5.246 

5.245 

5.124 

5.013 

5.245 

5.481 

5.509 

5.554 

5.097 

log(L/L©) 

1.909 

1.908 

1.588 

1.272 

3.376 

4.324 

4.439 

4.620 

1.286 

MHe/M© 

0.92281 

0.92249 

0.92205 

0.92204 

0.92198 

0.92156 

0.91634 

0.91600 

0.91600 

PHe 

11.594 

11.659 

4.156 

3.702 

0.724 

0.164 

0.171 

0.300 

7.417 

Thb 

1.248 

1.354 

4.073 

3.765 

3.949 

2.930 

29098 

2.812 

1.062 

log(LHe/L 0 ) 

3.052 

3.910 

11.311 

11.177 

8.767 

5.462 

5.319 

3.910 

-4.224 

A (in lO^^R©) 

9.651 

9.651 

9.651 

9.651 

9.651 

9.713 

9.718 

9.767 

9.615 

Macc/M© 

0.92985 

0.92985 

0.92985 

0.92985 

0.92985 

0.92307 

0.92253 

0.91707 

0.91707 

Mdon/M© 

0.13977 

0.13977 

0.13977 

0.13977 

0.13977 

0.13977 

0.13977 

0.13977 

0.13977 

At (yr) 

24468.8 

11.4436 

3.37408 

0.00006 

0.00325 

0.83678 

3.34515 

45.7288 

6028.94 


AMtran = 10.230 

AMio,t = 12.787 

Pacc = 

=-0.250 






SI 02+020 






log(Tefi) 

5.479 

5.479 

5.368 

5.375 

5.540 

5.607 

5.585 

5.629 

5.495 

log(L/L 0 ) 

2.744 

2.743 

2.477 

2.493 

4.361 

3.954 

4.541 

4.721 

2.808 

Mne/M© 

1.02090 

1.02078 

1.02068 

1.02070 

1.02051 

1.02070 

1.02010 

1.02036 

1.01983 

PHe 

4.560 

4.535 

1.395 

1.445 

0.455 

0.716 

0.468 

0.133 

2.754 

Thb 

1.488 

1.637 

4.130 

4.081 

3.543 

3.976 

3.249 

2.599 

1.612 

log{LHe/L©) 

3.867 

4.744 

9.573 

9.566 

7.175 

8.588 

5.178 

4.396 

-0.088 

A (in 10“^R©) 

8.020 

8.020 

8.020 

8.020 

8.020 

8.020 

8.029 

8.030 

8.016 

Macc/M© 

1.02250 

1.02250 

1.02250 

1.02250 

1.02250 

1.02250 

1.02100 

1.02067 

1.02067 

Mdon/M© 

0.19797 

0.19797 

0.19797 

0.19797 

0.19797 

0.19797 

0.19797 

0.19797 

0.19797 

At (yr) 

939.911 

0.81022 

0.26848 

- 0.00001 

0.05024 

-0.04771 

2.90204 

13.4218 

181.015 


AM 

tran = 2.034 

AMiost = 1.829 

Pace- 

= 0.101 






S060+017 Cold 





log(Teft) 

5.128 

5.127 

4.876 

4.626 

4.573 

4.388 

4.972 

5.440 

4.850 

log(L/L 0 ) 

1.883 

1.880 

1.195 

0.450 

0.534 

0.660 

3.399 

4.027 

0.642 

Mhb/M© 

0.60318 

0.60230 

0.60089 

0.60088 

0.60087 

0.60061 

0.59746 

0.59540 

0.59050 

PHe 

7.907 

7.883 

2.154 

1.258 

0.757 

0.186 

0.169 

0.266 

4.113 

Thb 

1.233 

1.341 

3.271 

3.558 

3.596 

3.022 

2.304 

1.664 

0.568 

log{LHe/L©) 

2.964 

3.882 

10.317 

9.770 

9.204 

6.714 

5.605 

2.716 

-7.129 

A (in 10“^R©) 

8.626 

8.627 

8.627 

%.(ai 

8.627 

8.667 

8.812 

8.958 

8.520 

Macc/M© 

0.62646 

0.62647 

0.62647 


0.62647 

0.62297 

0.61032 

0.59803 

0.59803 

Mdon/M© 

0.14032 

0.14031 

0.14031 

0.14031 

0.14031 

0.14031 

0.14031 

0.14031 

0.14031 

At (yr) 

58331.2 

45.0383 

11.9684 

0.0005 

0.0019 

0.2854 

30.3768 

66.2967 

86893.4 


AMt 

Iran — 29.694 

AMiost = 28.444 

Pace- 

=0.042 
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the accretors occurs inside “the Roche lohe finite space”, so that all 
the matter passing through the critical Roche surface is lost by the 
stars and, hence, by the binary systems. At the onset of the RLOF 
episode, due to the loss of matter and angular momentum from the 
systems, the orbital separation increases, therefore, the donors de¬ 
tach from their Roche lobes and, hence, the mass transfer halts. 
During the RLOF episode the accretors evolve at almost constant 
radii, increasing their surface luminosity and effective temperature 
up to when they recede definitively from their Roche lobe (point 
Ft). Since the components of the binary are now detached, so that 
mass transfer from donors to accretors can resume only after an¬ 
gular momentum loss by GWR shrinks the orbits, thus forcing the 
donors to overfill once again their Roche lobe. During the time be¬ 
tween the end of the RLOF and the re-onset of mass accretion, 
accretors evolve first up to the bluest point along the loop in the 
FIR diagram and, then, down along the cooling sequence. 

For each considered system, in Table we report the total 
mass transferred from the donor to the accretor during the first mass 
transfer episode (AMtran), the mass lost during the RLOF phase 
(AMiost) and the corresponding accumulation efficiency r^acc, de¬ 
fined as in Eq. (5) in Paper I: 

AMiost 

“ AMt.i + AMt.2 ’ ^ ' 

where AMtri is the amount of mass transferred up to the onset of 
the RLOF episode and AMtr 2 the one accreted after the end of the 
RLOF up to the bluest point along the loop in the FIR diagram. 


3.1 S060-I-017 System 

In Fig. m we plot for the S060-I-017 and S092-I-015 systems (the lat¬ 
ter to be discussed in the next subsection) the dependence of the 
mass transfer rate from the donor (upper panel) and the depen¬ 
dence of the slope of the curves n = d log M/d log Pori) (lower 
panel) on the orbital period Porb- For the sake of comparison, in 
the lower panel we show also the n curve for a system with initial 
parameters equal to the ones of S060-I-017, but evolving completely 
conservatively. The Figure clearly shows that thermonuclear out¬ 
bursts, interrupting the mass transfer (see previous Section), virtu¬ 
ally do not influence the common behaviour of the M- Porb rela¬ 
tion for ultracompact binaries evolving un der the influence of an¬ 
gular momentum loss by GWR. Recently, ICannizzo & Nelemanj 
( 1201 SI) , using the limit cycle accretion disk instability model, at¬ 
tempted to constrain semi-analytically the mass transfer rate in 
AM CVn-stars by fitting M to the Porb limits of the accretion 
disk-instability range of these systems, namely 20 and 44 min. 
The metallicity of AM CVn-stars is unkn own. For the range of Z 
from 0.0 to 0.04, ICannizzo & Neleman j found values of n from 
—5.38 ± 0.3 to —5.06 ± 0.3. Though the period range of interest 
fo r the present work is below 20 min., the comparison to the results 
of ICannizzo & NelemansI (12015 ), as well as to earlier evalua tions 
of, e. g., Tutukov & YungelsonI 1996h : l^lemans et~^ ( l200lh val¬ 


idates our computations. As it concerns the estimates of accretion 
rates in AM CVn-stars, they are obtained by indirect methods and 
hardly can constrain the M-Porb relation. They may be considered 
as “indicators” of the proper order of magnitude of the computed 
M, at best. 

In Fig.|5]we plot the mass transfer rate (upper panel) for the 
S060-(-017 system as well as the evolution of the temperature (mid¬ 
dle panel) and density (lower panel) of the He-shell in the accre¬ 
tor as a function of time. As it can be seen, after the He-flash and 
the related RLOF episode, mass transfer is not active and the He- 



4 6 8 10 20 

Porb (min) 

Figure 4. Upper panel: dependence on the orbital period of the mass trans¬ 
fer rate from the donors in the systems S060-(-017 (solid line) and S092+015 
(dashed line). The small vertical ticks on the curves mark the values of 
Porb when mass transfer does not occur. In these systems the first out¬ 
bursts happen at Por()=4.792 min. and 4.830 min., respectively and the cor¬ 
responding gaps APorb (0.13 and 0.09 s, respectively) are too short to be 
noticeable in the scale of the plot. Filled circles with errorbars show es¬ 
timates of mass transfer rates in obs e rved A M CVn systems (taken from 
the the compilation by iKotko et*^ j2012h i. Lower panel: the value of 
n = d log M/d log for the S060+017 and S092+015 systems (solid 
and dashed lines, respectively). The dotted line shows the n profile for the 
case of completely conservative mass exchange in a (0.60+0.17) Mq sys¬ 
tem. 


burning shell cools down. When mass transfer resumes, the cor¬ 
responding accretion rate is about 3.23 x 10~^ M© yr“^. Due to 
the deposition of matter, the He-shell starts to heat up once again, 
but when M becomes lower than 8 x 10“® M© yr“^, the mass 
deposition becomes unable to balance the cooling and the accretor 
further evolves along its cooling sequence. Thus, despite quite sub¬ 
stantial amount of matter may still be transferred to the CO WD in 
this system (about O.ITM©), the nuclear evolution of the binary is 
terminated. Though, it will retain its status of AM CVn star. 

At variance with the results discussed above, the “heated” 
M060 model in Paper I, corresponding to the accretor in the 
S060+017 system, after the first He-flash, burns helium quiescently 
for M^ 1.5 X 10“^ M© yr“’^. However, it has been noted that 
the present computation differs from those in Paper I in many re¬ 
gards. First of all, in Paper I the accretion rate was kept strictly 
constant in each simulation and mass deposition was restarted soon 
after the end of the RLOF, thus preventing substantial cooling of 
the He-buffer. In addition, the adopted value for the Roche lobe 
radius (10 R©) was definitively larger than the one proper to the 
self-consistent computation of the S060+017 system. 
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Figure 5. Time evolution of the accretion rate (M - upper panel), of the 
He-shell temperature (Tjje - middle panel) and density (pHe ■ lower panel) 
of the accretor in the binary system S060+017. Dotted lines refer to the 
evolution of the same object, but starting to accrete mass soon after the end 
of the RLOF. For more details see text. 

In order to illustrate the origin of such a difference, we com¬ 
puted two models, arhitrarily varying the value of the separation 
during the evolution of the S060+017 system. In the first one we 
reduced the value of the separation after the end of the RLOF 
episode in such a way that mass transfer resumed immediately. 
The time evolution of the mass transfer rate and of the tempera¬ 
ture and density of the He-burning shell for this model is displayed 
by dotted lines in Fig.O The prompt onset of mass accretion pre¬ 
vents the decrease of the temperature of the He-buming shell below 
~ 8.3 X 10^ K, but, in any case, as M continuously decreases, the 
accretor never attains the physical conditions suitable for the re¬ 
ignition of helium. This numerical experiment clearly suggests that 
the origin of the difference between the present computation and 
the corresponding one in Paper I depends on the different thermal 
content of the post-RLOF structure. 

To make this more clear, we performed the second run (Test 
model) increasing arbitrarily by a factor of 100 the binary sep¬ 
aration at the epoch of the onset of the RLOF overflow in the 
S060+017 system. Note that, as a consequence of such a variation, 
the Roche lobe of the accretor became ~ 4R©, while the Roche 
lobe radius of the donor increased to 2.2 R© and, hence, the mass 
transfer came to a halt. In this way we constructed a model having 
exactly the same mass of the He-rich buffer of the accretor as in the 
S060+017 system, but which “freely” expanded in the space before 
the onset of the RLOF, as in the computations of Paper I. 

In the upper panel of Fig.l^we show the time evolution of the 
flash-driven convective shell (gray area), of the total mass (dotted 



-50 0 50 


Age-Agenu,p [yr] 

Figure 6. Upper panel: time evolution of the flash driven convective shell 
(gray zone), of the total mass (dashed line) and of the mass coordinate of the 
He-buming shell (solid line) for the accretor in the 8060+017 system. The 
vertical dashed line in the upper panel marks the end of the RLOF episode. 
Lower panel: behaviour of the same quantities at the same time in the Test 
model where RLOF occurs later. 

line) and of mass coordinate of the He-buming shell (solid line) 
for the accretor in S060+017 system. In the lower panel we show 
the behaviour of the same quantities, but for the Test model, which 
has larger Roche lobe. The epoch t=0 corresponds to the onset of 
the RLOF episode in the S060+017 system, while the dashed ver¬ 
tical line marks its end. As it is seen, the mass loss triggered by the 
RLOF episode determines the rapid disappearance of the convec¬ 
tive shell. As discussed in Paper I, convection affects the evolution 
of the He-flash in two opposite ways: on one hand it removes en¬ 
ergy from the inner zones of the He-rich buffer and redistributes it 
over the whole envelope, thus reducing the local increase of tem¬ 
perature and damping the thermonuclear runaway. On the other 
hand, convective mixing dredges down fresh helium, feeding the 
burning-shell and, hence, powering the flash. As a result, the strong 
and rapid reduction of the convective shell in the S060+017 sys¬ 
tem limits the amount of nuclear energy produced via He-buming 
and, hence, the resulting flash is weaker with respect to the Test 
model. To make this conclusion more quantitative, we define for 
the model S060+017 two time intervals: Afi = 55.468 yr, lasting 
from the onset of the flash-driven convective shell to the onset of the 
RLOF, and At 2 = 69.774 yr, lasting from the onset to the end of 
the RLOF. In the S060+0I7 model during Ati He-buming delivers 
an amount of energy equal to £i ~ 2.678 x 10^® erg, while during 
At 2 it delivers £2 — 1.639 x 10^® erg. In the Test model, the energy 
delivered during At 2 is by a factor 2 larger (£ 2 ^^ — 3.219 x 10"*® 
ergfl 

Figure 1^ also demonstrates that the extension and lifetime of 
the flash-driven convective shell affects the heating of the layers 
below the He-buming shell. This is clearly indicated by the inward 

Note that during At 2 , the surface radius of the accretor in the 8060+017 
system remains practically unaltered, as it is fixed by the Roche lobe radius, 
while in the Test model it increases by a factor 25. 
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Figure 7. Dependence of the total luminosity L (upper panel) and ratio of 
luminosity of the He-buming shell Ljje and total luminosity L (lower panel) 
on the radius of the CO WD in the sequence S062+017 during the He-flash. 
S062+017 system is shown with solid lines, Test model with dashed ones. 
In the evolutionary phase where the two accretors evolve identically, dashed 
lines are hidden behind the solid ones. Different points and letters mark 
various important epochs, as in Fig.|2] 


Figure 8. Comparison of the properties of the first flash in the S060+017 
(dashed lines), S092+015 (solid lines), and S102+020 (dotted lines) sys¬ 
tems. The panels (top to bottom) show the luminosity of the He-buming 
shell of accreting dwarf, the total energy released during the flash, and the 
“specific energy” of the flash (see text for details). The age of the S102-1-020 
system has been arbitrarily increased by 1.9 X 10^ yr. 


shift of the He-buming shell (see solid lines in Fig. | 6 j: a larger 
amount of thermal energy produced via He-burning determines a 
more efficient heating of the underlying zones and, hence, a deeper 
inward shift of the He-burning shell itself. 

To illustrate further the difference between these two cases, i. 
e., between the evolution of very compact systems and that of wider 
ones, we present in Fig.|7]the dependence of the total luminosity of 
the accretor L (upper panel) and the ratio of the He-buming shell 
luminosity and the total luminosity Lne/L (lower panel) on the 
WD radius. The evolutionary paths of two accretors coincide up 
to RLOF epoch in the S060-I-017 system. In the latter, soon after 
the onset of the RLOF, convection is aborted and the He-rich lay¬ 
ers above the He-buming shell, previously heated by the flash, are 
lost in a short time scale (Fig. [D- This almost extinguishes He- 
burning very rapidly and the only energy source to maintain ther¬ 
mal equilibrium becomes contraction. This can be seen in the lower 
panel of Fig.|7] showing that after the RLOF episode in the system 
S060-I-017 the ratio Lne/L becomes ~l/20. As He-buming dies 
completely, contraction will remain the only energy source balanc¬ 
ing the radiative losses from the surface since He is never reignited 
because of the cooling of the He buffer due to the interruption of 
mass-transfer and, later on, to the insufficient heating via accre¬ 
tion (Fig. |5j- In the Test model, convection continues to support 
the He-burning shell and, consequently, the delivered energy drives 
the expansion of the WD He-rich envelope to the large Roche ra¬ 
dius. The dwarf recedes from its Roche lobe when the energy pro¬ 
duced via He-buming equates the energy losses from the surface. 


However, this occurs slower than in the case of the tight system - 
for the same luminosity level after the RLOF, the accretor in the 
S062-I-017 model is more compact then in the Test model (upper 
panel of Fig. |7j. Later on, the two models converge to the same 
evolutionary path. 

The difference illustrated above has a sizable effect on the re¬ 
tention efficiency of the two models. In particular, in the Test model 
a larger amount of nuclear energy is delivered before the onset 
of the RLOF episode, so that the the whole He-buffer above the 
CO core is heated more than in the S060-I-017 system. As a con¬ 
sequence, the He-burning shell can move inward; moreover, even 
if a part of the energy injected in the He-rich buffer is dissipated 
via mechanical work (expansion to larger radius), the thermal con¬ 
tent of the He-rich mantle remains very large. As a matter of fact, 
a larger amount of matter is lost during the RLOF episode in the 
Test model (AMiost = 2.7 x 10“®Mq), so that the corresponding 
retention efficiency reduces to rja,cc=0.065 (about 45% lower than 
the one obtained for the S060-I-017 system). In Paper I we found 
that the “heating” flash in the M060 Cool model erodes also part 
of the pre-existing He-rich buffer, while in the Test model we find 
a small, but in any case positive, retention efficiency. Such a dif¬ 
ference reflects the different values of M adopted for the above 
mentioned models: in Paper I we use M=10“^ Mq yr“^, while in 
the Test model the accretion rate is definitively larger: initially it is 
as high as 8.8 x 10 “^ M© yr“^ and, in any case, it remains larger 
than 3.3 x 10“^ Mq yr“^). In Paper I we explored the dependence 
of r)acc on the assumed Roche lobe radius, by varying the latter in 
the range 1-45 R 0 and we found that the difference in the estimate 
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Table 3. Selected physical properties for each of the five strong pulses ex¬ 
perienced hy the accretor in the S092-I-015 system. Mhb (the mass coordi¬ 
nate of the He-buming shell, Mq ), AMjje (the mass above the He-buming 
shell, 1 O“®M 0 ), Porb (the orbital period, in min), pjje and Thb (density 
in lO^g • cm“® and temperature in 10® K at the He-burning shell) refer 
to the epoch of He-ignition. is the maximum luminosity of the He- 

burning attained during the He-fiash. is the mean value of the mass 

accretion rate (in 10“^ Mq yr“^) during the time lasting from the onset 
of mass transfer and the onset of the RLOF episode. AMtran and AMio^t 
(in 10“®Mq) are the amounts of mass transferred from the donor and lost 
by the accretor during the RLOF, respectively, rjncc represents the accumu¬ 
lation efficiency, while At (in yr) is the time lasting from the end of the 
RLOF to the re-onset of mass transfer. 




2Ti.fi 

yd 

yh 

yh 

Mhc 

0.9225 

0.9180 

0.9178 

0.9189 

0.9203 

AMhc 

7.362 

5.887 

6.851 

8.935 

14.170 

Porb 

4.830 

5.069 

5.331 

5.753 

6.588 

PHe 

11.66 

9.22 

10.78 

14.00 

22.13 

Thb 

1.354 

1.385 

1.367 

1.289 

1.205 

log(L“-VL0) 

11.311 

10.787 

10.967 

11.592 

12.520 

(m) 

4.178 

3.087 

2.413 

1.748 

1.015 

AMtran 

10.230 

6.865 

6.838 

9.799 

15.802 

AMiost 

12.787 

6.107 

6.584 

9.191 

15.085 

^acc 

-0.250 

0.110 

0.037 

0.062 

0.045 

At 

6028.9 

3434.5 

4805.2 

8401.0 

23316.8 


of retention efficiency is smaller than 7-8%. Such a conclusion is 
appropriate for binaries with relatively massive He-star donors, but, 
as demonstrated above, not for low-mass compact stars. 

To summarize, in the S060-I-017 system, due to the small sepa¬ 
ration, the first He-flash does not result in an efficient heating of the 
layers above the CO core, so that the following evolution proceeds 
substantially different from what derived in Paper I. 


3.2 S092-I-015 System 

The first flash in the S092-I-015 system is much stronger than in 
the other systems, as clearly displayed in Fig.[8l where we plot as 
a function of the evolutionary time three quantities characterising 
the He-flash for all the systems listed in Table the He-burning 
luminosity Lhb (upper panel), the total energy released during the 
flash Ehs (middle panel), and the “specific energy” of the flash eee 
(lower panel), defined as the ratio of Ehb and mass of the He-rich 
zone. 

The maximum luminosity in the burning shell Lhb attained 
during the flash is 10 times larger than for the S060-b017 system 
and about 50 times larger than for the SI02-1-020 one (see also Ta¬ 
ble |3. Such an occurrence is not related to the physical properties 
of the He-shell and to the thermal content of the underlying CO core 
at the onset of the mass transfer, but is determined by the adopted 
combination (Mdon, Mace) which sets the value of M after contact. 
At the beginning of the mass transfer, M for the S092-I-015 system 
is ~ 5.0 X 10“^ while for the S060-I-017 and S102-1-020 systems 
it is 8.8 X 10“^ and 2.2 x 10“®, respectively. As a consequence, 
the compressional heating of the He-buffer occurs at a lower rate 
so that at the onset of the He-flash the He-shell has more degen¬ 
erate physical conditions (see the values of pHe and Thc listed in 
column B of Table|^. 

After the first He-flash, at the re-onset of mass transfer M~ 
3.5 X 10“^ Mq yr“^; according to Paper 1, the same initial mass 
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Figure 9. The same as in Fig.[5]but for the S092+015 system. 
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Figure 10. Comparison of the properties of subsequent flashes in the se¬ 
quence S092-I-015. The same variables as in Fig.[ 8 ]are plotted. 
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WD, model M092, after the “heating flash” experiences strong He- 
flashes for 5 x 10“^ Mq yr“^. The S092+015 system un¬ 
dergoes other 4 strong non-dynamical He-flashes before a massive 
degenerate object is formed, (see Table 0 and Figs. |9TT0 ] i. In Ta¬ 
ble [3] we list selected physical quantities of the accretor during the 
5 He-flashes experienced by this system. 

By adopting for M the mean value before the flash re¬ 
ported in Table [3 and interpolating the data reported in Table 4 of 
Paper I we obtain ?7acc~ 0.61, definitively larger than the value 
we derive in the present worl|3- As discussed in the case of the 
S060-I-017 system, such a discrepancy has to be ascribed to the 
non-efficient heating of the He-layer during the first He-flash due 
to the sudden and sharp decrease of the mass extension of the flash- 
driven convective shell. When mass transfer resumes after the 5*^ 
flash, M= 7.25 x 10~® M© yr“^ and the He-shell starts to heat 
up once again. At M~ 4 x 10“® Mq yr“^ the accretor enters the 
regime which for the actual Mwd and constant M would corre¬ 
spond to dynamical flashes (Paper I). In the case under analysis, 
as M continuously decreases, neutrino cooling starts to dominate 
over the compressional heating and no additional He-flashes (dy¬ 
namical or not) are ignited. Hence, accretor turns into a nuclearly 
inert degenerate object which gradually increases its mass. 

3.3 S102-I-020 System 

The evolution of the S102-1-020 system after the first flash is quite 
different. Figure and Table reveal that for this system the evo¬ 
lution during the first He-flash episode is more similar to that of 
a “freely” expanding model. In fact the flash driven convective 
shell attains the stellar surface before the epoch of the maximum 
He-shell burning luminosity and, in addition and more important, 
the RLOF episode starts after the flash-driven convective zone has 
receded from the stellar surface (see the sequence of the various 
phases in Fig. [3 and the negative value of Af for epochs D and F 
in Table O. This means that the He-flash succeeds in heating the 
whole He-buffer to a higher level with respect to what occurs in the 
S060-I-017 and S092-I-015 systems. Such an occurrence is a direct 
consequence of the fact that, when the components in this binary 
system come to contact for the first time, the mass transfer rate is 
very high (M~ 2.2 x 10“® Mq yr“^), so that gravitational en¬ 
ergy released by accretion efficiently heats up the entire He-buffer. 
This determines less degenerate physical conditions at the onset of 
the He-flash, as compared to the other two considered systems (see 
PHe and Ths for the epochs A and B in Table[2j, so that the result¬ 
ing He-flash is less strong (see the value of Lhb at the epoch C in 
Tableland Fig.[lll. 

When mass transfer from the donor resumes, M is still as high 
as 2.1 X 10“® Mq yr“^(see upper panel in Fig. [IB. In Pa¬ 
per I we found that the “heated” MI02 model accreting at such 
a rate experiences quiescent burning of He. On the contrary, in 
the present work the accretor in the S102-1-020 system experiences 
other 3 flashes. The flashes become progressively less strong (see 
middle panel of Fig. [IB and the corresponding retention efficiency 
increases (see lower panel of the Fig. [IB- In particular, we find 


® According to the discussion in Sect. 3.1 of Paper I, the evolutionary out¬ 
come of WDs accreting mass with a time-dependent accretion rate depends 
mainly on the current value of M and to a less extent on the previous ther¬ 
mal history. At the onset of the 2"'* He-flash in the S092-t015 system we 
find Mai 2.5 X 10“^ Mq yr“^; for this M, according to the data in Table 
4 of Paper I, we obtain 77acc—0.49. 



-3 - 2.5 -2 - 1.5 


log(0.2-M,„„/M3) 

Figure 11. Evolution of the S102-t020 system. We report as a function of 
the mass transferred from the donor the mass transfer rate M (upper panel), 
the luminosity of the He-buming shell (middle panel) and the mass of the 
accretor Mace (lower panel). Vertical dashed line in the upper panel marks 
the transition from one accretion regime to another (see text for details). 


?yacc=0.562 for the second flash and r7acc=1.0 for the third one. We 
classify the latter as a “mild flash” (MF), i. e., a flash which re¬ 
leases so little nuclear energy that the accretor remains confined 
well inside its Roche lobe and no mass-loss occurs. Such a be¬ 
haviour is determined by the fact that, pulse by pulse, the He-rich 
buffer and the most external layers of the underlying CO core heat 
up, thus attaining the physical conditions suitable for quiescent He- 
buming. The resulting evolution is completely different from that 
of the accretor in the S092-I-015 system and it is a direct conse¬ 
quence of the mass transfer rate after the first flash episode which 
depends on the parameters of the initial binaries. In particular, as 
already mentioned before, M in the S092-I-015 system after the 
first flash should determine recurrent strong flashes also in a fully 
heated model with the same total mass and mass of the He-huffer 
(e. g. the “heated” M092 model in Paper I). At variance, in the 
SI02-1-020 system the mass transfer rate is so high that the released 
gravitational energy prevents the cooling down of the He-buming 
shell during each inter-flash period. As displayed in Fig. [IB af¬ 
ter four flashes the accretor reaches Steady burning regime (Steady 
State), which lasts as long as, due to continuous decrease of the 
mass transfer rate from the donor, the extension of the He-buffer 
above the He-burning shell reduces below a critical value and the 
accretor enters again the mild flashes regime. The transition from 
one regime to another occurs smoothly so that we arbitrarily de¬ 
fine the value of M at which such a transition occurs as the epoch 
when the maximum luminosity of the He-buming shell during the 
He-flash becomes twice the value of the surface luminosity along 
the high luminosity branch. Under this assumption we find M(SS- 
MF)~ 1.09 X 10“® Mq yr“^. Due to the evolution of binary pa¬ 
rameters, the mass transfer rate continuously decreases and when it 
becomes lower than ~ 8.3 x 10“^ Mq yr“^, after 29 MFs, Strong 
Flashes start again. In the upper panel of Fig. m the transitions 
from one regime to another are marked by dashed vertical lines. 

The following evolution is described in Fig. [IB where we 
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Figure 12. Evolution of the accretor in the Si02+020 system during the 
Strong Flashes accretion regime. We report, as a function of the total 
mass of the accretor at the bluest point along the loop in the HR diagram 
(Myyu blue) Value of the retention efficiency (panel a), the luminosity 
of the He-buming shell (panel b) and the temperature (panel c) and density 
(panel d) at the epoch of He-ignition. 

show, as a function of the accretor total mass at the bluest point 
along the loop, the value of the retention efficiency, and the physical 
conditions (temperature and density) at the epoch of the He-flashes 
ignition. As well, the maximum luminosity in the He-burning 
shell Lhb is shown. As expected for a monotonically decreasing 
mass transfer rate, ?7acc, on average, decreases continuously (upper 
panel), as the physical conditions at the He-shell become more de¬ 
generate and, hence, the strength of successive flashes increases. 
The general trend is defined by two main factors. On one hand, 
as the accretion rate decreases and T^acc reduces, the time span af¬ 
ter the end of the RLOF episode and the mass transfer resump¬ 
tion increases, so that the the He-buffer cools more. This deter¬ 
mines the increase of the amount of mass to be transferred to ignite 
the successive He-flash which, as a result, is also stronger. On the 
other hand, pulse by pulse, the Roche lobe radius of the accretor 
increases, so that during the successive He-flashes the heating of 
the whole He-buffer is more efficient, as the RLOF occurs when 
a larger amount of nuclear energy has been delivered (see above 
the discussion of S060+017 system and related Fig. 1^. Though, 
the evolutionary curves also exhibit irregularities whose origin de¬ 
pends on the complex interplay of various additional factors. 

We find that after the 58*^ strong flash, the mass of the donor 
has reduced to ~O.O875M0 and that of the accretor increased to 


~ 1.0822Mq. When accretion resumes, after 181916 yr, the mass 
transfer rate is 4.2 x 10“® Mq yr“^ and it rapidly decreases so that 
the compressional heating produced by accretion is not able to bal¬ 
ance the radiative and neutrino cooling of the whole He-rich buffer. 
Hence, as already illustrated for the S060+017 and S092+015 sys¬ 
tems, the accretor in the S102+020 system cools down and the fi¬ 
nal outcome is a massive CO core surrounded by an extended He- 
buffer. 


3.4 “Cold” S060+017 Model 

In order to investigate the dependence of our results on the cool¬ 
ing age of the accretor, we let the “heated model” M060 from 
Paper I to cool for 2 Gyr which, according to population synthe¬ 
sis computations, can be considered as the typical time of for- 
mation of an AM CVn sy stem from a pair of detached WDs 
jTutukov & Yungelsoiilll99^ . After that, we start the accretion fol¬ 
lowing the same procedure as above. The results are summarised in 
Table 12 As it can be noticed, at the beginning of the mass transfer, 
the physical conditions at the base of the He-burning shell are more 
degenerate (temperature is a factor 5 lower, while density is practi¬ 
cally the same). Notwithstanding, the amount of mass to be trans¬ 
ferred to ignite the He-flash is only 3% larger and the ignition con¬ 
ditions are very similar (compare columns A and B in Tabs. |2 for 
the standard and “Cold” S060+017 Models) and the resulting He- 
flash has almost the same strength. During the RLOF episode the 
cold S060+017 system loses about 11% more mass so that the fi¬ 
nal retention efficiency is lower. As a consequence, the post-RLOF 
episode separation is a bit larger and, hence, the time span up to the 
re-onset of mass transfer increases by more than a factor 3.5. In any 
case, the evolution after the He-flash episode occurs exactly as in 
the S060+020 system: the mass transfer rate decreases very rapidly, 
so that the radiative and neutrino losses from the He-buffer become 
dominant and the system evolves to the formation of a massive CO 
WD with an extended He-rich mantle. 

The results for the “Cold” model suggest that the models we 
considered are representative also for systems with longer cooling 
age. The reason is that at the luminosity level 0.01 Lq the physical 
conditions at the center of WD and at the base of the He-rich buffer 
have attained an asymptotic value. Increasing the cooling age of ac¬ 
creting WDs from 800 Myr to (2 - 3) Gyr practically does not affect 
their reaction to accretiorQ. Such a statement is also confirmed by 
the computation of the evolution of the “Cold” analog of S092+015 
system up to the onset of the first RLOF episode. The released nu¬ 
clear energy was quite similar to the one in the “hot” model, while 
the accreted mass was only 1.05% larger. This similarity suggested 
us to terminate this additional extremely time-consuming computa¬ 
tion, since we infer that the initial temperature of WD will not play 
any role, like in the case S060+017. 


4 DISCUSSION AND CONCLUSIONS 

In the present work we studied accretion from He WD donors onto 
CO WDs in ultracompact AM CVn binaries (IDD). At variance 
with our earlier study of the He-buming regimes in WDs accreting 
mass at constant M (Paper I), in our current analysis we adopted 


° Longer formation times are not considered in the present analysis as they 
are not typical for AM CVn stars. 
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time-dependent accretion rates, as determined by the loss of angu¬ 
lar momentum via GWR. As in Paper I, we assumed that, if He- 
flashes on the accretors result in RLOF, matter is lost from the sys¬ 
tem until accretor recedes from the critical lobe. The matter leav¬ 
ing the system has the specific orbital angular momentum of the 
accretor. The systems considered in our analysis have the follow¬ 
ing properties: (i) mass transfer rate permanently declines and (ii) 
accretion may be interrupted due to mass and angular momentum 
losses. In Paper I we assumed that the accretion process is almost 
continuous, as it resumes soon after the end of the RLOF episode. 
At the end, in our earlier study we preset arbitrarily the accretor 
Roche lobe radius, while in the present work the latter is defined 
by the actual parameters of the given binary system, i. e., masses of 
the components. The initial values of them were taken correspond¬ 
ing to supposed precursors of AM CVn stars — extremely low 
white dwarfs systems (ELM WDs). The donors were approximated 
as pure He-objects with zero-temperature, proper to the formation 
timescale of AM CVn stars ~Gyr. 

The actual He-burning regime onto the accretors in IDD is 
defined by the interplay of several factors. The first, as in other ac¬ 
creting systems, is the balance between the compressional heating 
driven by the mass deposition and the cooling of the He-rich buffer 
via inward thermal diffusion and neutrino emission. The second is 
the permanently decreasing accretion rate. The third is the degen¬ 
eracy level of the physical base of He-rich buffer. Last but not least, 
the Roche lobe radius of accretors may be very small, down to sev¬ 
eral 0.01 R©. 

As expected (see the discussion in Paper I), we found that af¬ 
ter the onset of mass transfer onto the accretors, the He-burning 
shell above the CO core heats up and, when the nuclear timescale 
for Sa-reactions becomes shorter than the local thermal adjustment 
timescale, a nuclear flash occurs. Due to the degenerate physical 
conditions, this first flash is strong and it drives to a significant ex¬ 
pansion of the star. It is worth noticing that this first flash is the 
equivalent in the real world of the “artificial” initial heating flash 
used in the studies of He-buming onto WDs, in order to mimic post- 
AGB objects. However, in a typical AM CVn system the accretor 
should be cold, since the timescale of AM CVn stars formation is 
~ Gyr. As a matter of fact, such a flash results in a RLOF by the 
accretor and drives to the loss of mass and angular momentum from 
the system. 

As shown in Sect. 13.11 at variance to “wide” interacting sys¬ 
tems, the post-flash evolution of IDD is affected by the prompt 
mass loss which limits the lifetime and extension of the flash-driven 
convective shell, thus limiting the feeding of the He-burning shell 
by convection and reducing the total amount of nuclear energy re¬ 
leased during the flash. Such an occurrence prevents the efficient 
heating of the pre-existing He-rich layers below the He-burning 
shell and causes the rapid extinction of nuclear burning. Hence, 
in the post-RLOF phase the accretors in IDD contract very rapidly 
in order to maintain the thermal equilibrium. 

According to the accepted paradigm for the evolution of out- 
bursting binaries, the matter overflowing the accretor Roche lobe 
leaves the system, taking away the specific angular momentum of 
the accretor. As a result, the system becomes detached, mass trans¬ 
fer comes to a halt and the WD cools. Hence, the possibility of 
resumption of nuclear activity depends on the duration of this “hi¬ 
bernation” phase and the masses of the two components which, in 
turn, determines the value of M at the the resumption of mass trans¬ 
fer. As it has been found for the sequence S062-I-017 (Sect. 13. il l, for 
an ever decreasing M, it is possible that the release of gravitational 
energy occurs in a so long timescale that it can not counterbalance 


the cooling by inward thermal diffusion and neutrino emission and 
another flash can not be ignited. For the other systems considered in 
the present work, this limits the total number of flashes experienced 
by a given system. 

The rough comparison of the present study and Paper I sug¬ 
gests that the expected number of nuclear flashes in IDD is a factor 
2 lower than the one expected by interpolating the results based 
on models with constant M. This reduces further the estimation ob¬ 
tained in Paper I of only several thousand “nuclearly active” IDD 
currently present in the Galaxy. 

The results obtained in the current work exhibit a very mild 
dependence on the initial thermal content of the accretor, for cool¬ 
ing ages in the range 0.8 - 2 Gyr (see Sect. l3.4b . so that they can be 
considered as representative of the entire population of IDD. 

For each fixed initial CO WD mass, as the mass transfer rate 
is permanently decreasing, the last flash experienced by accreting 
WD s is the strongest one (see Figs.l^ [70lll2b . As mentioned be¬ 
fore, [siidstenlet^] l l2007h suggested that in a series of flashes the 
last one may be of dynamical nature (SN .la). None of the sys¬ 
tems considered in the present work develops the physical condi¬ 
tions suitable for such a dynamical flash. In particular, we found 
that the mass of the He-buffer at the onset of the last He-flash is 
always lower than the critical value necessary to obtain a SN .la, 
derived by Bildsten and his c oauthors for the same WD mass (se e 
Fig. 2 in [Bildsten et al.|[200^ and Fig. 1 in ishen & Bildstenll2009h . 
Such a difference has to be ascribed to the different thermal con¬ 
tent of the layer below the He-burning shell at the onset of the 
last He-flash. In our model, the previous accretion history deter¬ 
mines the partial heating of the most external zone of the CO WD 
and of the He-buming shell, so that at the beginning of the last 
mass transfer episode driving to a He-flash we have Thb ~ 10® K, 
for both S092-I-015 and S102-I-020 systems. As a consequence the 
“last flash” in our computations will never attain physical condi¬ 
tions suitable to synthesize iron peak elements. Thus all three con¬ 
sidered systems have the same fate: transformation of the accretor 
into a massive degenerate object with a CO-core, which has almost 
the same mass as the initial accretor, and a massive He-buffer, rep¬ 
resenting a lar ge part of the i nitial d onor. 

Recently, [Brooks et alj (1201 Sh performed an analysis similar 
to our one, but focused on AM CVn systems with a non-degenerate 
He-buming star as a donor. Their results can not be compared di¬ 
rectly with our findings, because the response to mass extraction of 
a He-star is different from that of a degenerate object and, hence, 
the resulting evolution o f the host binary system has to be differ¬ 
ent. [Brooks et alj (1201 Sh found that the last He-flash experienced 
by the accretor of the “He-star family” of AM CVn systems is 
never strong enough to produce a SN .la event. So, by combining 
their results with our one, we can conclude that AM CVn family 
at large will never produce an explosive event at the end of their 
life. We remark also, that we follow in detail the thermal evolution 
of the accretors during each flash episode, determinin g the effec¬ 
tive accumulation efficiency, while [Brooks et al] ( 1201 Sh remove the 
entire He-envelope above the He-buming shell at the onset of the 
ther monuclear runaway. 

[Brooks et al.l ( l2015l) claim that the first He-flash experienced 
by AM CVn systems with non-degenerate donors more massive 
than 0.4M© and accretors more massive than 0.8 Mq is “vigorous 
enough to trigger a detonation in the helium layer”, which could 
produce either a SN .la or a real type la Supernova, if also the 
CO core detonates. Interestingly enough, this first explosive event 
could destroy the host binary, thus reducing the expected number 
of AM CVn systems from this evolutionary channel. This could 
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Figure 13. Mass transfer rate in the S092+015 system as a function of the 
mass transferred from the donor AMtran in different timescales t (in 
Myr), as labelled (solid lines). The dashed line represents accretion rate 
MpD corresponding to fully conservative mass transfer. Filled circle marks 
the coordinate {AMtran, M) at the onset of the thermonuclear runaway for 
the S092+015 system under “standard” assumptions. Filled squares repre¬ 
sent the same quantities for the models computed with different t. The two 
vertical lines mark the critical mass of the He-buffer for a 0.92Mq CO WD 
con-esponding to the c ondition Thmt = W Td,m. (left) and Theat = Tdyn 
(right) as derived from iBrooks et al.N2015h . 


help in solving the problem of the a pparent deficiency of observed 
AM CVn stars jCarter et alj|2013ah . At variance, we find that the 
first He-flash occurring in the IDD explored in the present study 
never develops the physical conditions suitable for a detonation in 
the accreted He-rich buffer. This difference reflects the different be¬ 
haviour of the mass transfer rate before the onset of the He-flash: in 
particular, for IDD M is contin uously decreasing in time, while it is 
rising in the binaries st udied by Brooks_et^ l l2015l) . In this regard, 
it is worth noticing that iBrooks et alj(l2015l) consider AM CVn sys¬ 
tems with He-stars unevolved or only slightly evolved at the onset 
of the RLOF. However, the mass transfer history from He-stars de¬ 
pends on the extent of the donor evolution prior to the RLOF. As a 
matter of fa ct the M curve ma y be rising, almost flat or very slightly 
decreasing iYungelsonll2Q08h . Moreover, material transferred from 
evolved He-star donors has lower He abundance reducing the prob¬ 
ability of strong outbursts. 

We recall that in our computation we model the He-WD donor 
as a zero-temperature object, i. e. as a completely degenerate one. 
In the real world, however, the most external layers of WDs, inde¬ 
pendently of their chemical composition, are only partially degen¬ 
erate or not degenerate at all. As a consequence, when they fill their 
Roche lobes, the mass transfer rate increases in time, as it occurs 
for homogeneous He-star. Once these surface layers have been lost. 


and the fully degenerate inn er region starts to be removed, M de¬ 
creases jPelove et alj|200'^ see also Fig. 13 in Paper I). The dura¬ 
tion of the first phase of mass transfer depends mainly on the mass 
of the non-degenerate layers and, in turn, on the cooling age of the 
donor. On a general ground, such an occurrence implies that the 
compressional heating of the accretor should occur, at least at the 
beginning, in a longer timescale, so that the resulting He-buffer at 
the onset of the He-flash should be larger. In order to evaluate quan¬ 
titatively the effect of this rising phase of the mass transfer rate, we 
computed several additional models for the S092-b015 system, by 
adopting as mass accretion rate: 

M = min • log (3) 

where Mp d is the rate derived with the M — R approximation de¬ 
scribed in Sect. 12 In the second term, the mass of the partially de¬ 
generate/not degenerate layers AMnd is parameterized by means 
of the timescale r: large values of r correspond to large AMnd- 
Mpp) — 4.978 X 10“^ Mq yr"*^ is the value of Mpp, at the on¬ 
set of the mass transfer. In Fig. [T2 we report as a function of the 
mass transferred from the donor the evolution of the mass trans¬ 
fer rate described by Eq. 13, for different values of r. This figure 
discloses that for realistic values of AMnd ^ 0.03 Mq, corre¬ 
sponding to r ^ 2 Myr, AMtran increases up to a factor of ~3, so 
that the resulting He-flash is stronger than the one in the S092-I-015. 
In any case, we do not expect that this could significantly alter the 
following evolution, because, as discussed in Sect. 13.11 the prompt 
occurrence of the RLOF limits the heating of the He-buffer and, 
hence, the effects on the thermal properties of the accretor. Values 
of AMnd ^O.OSMq, corresponding to r > 5 Myr, are rather un¬ 
likely. According to the results of these numerical experiments, we 
can conclude that a dynamical event could hardly arise in IDD. 

Helium WDs with masses in the range explored in the present 
study have non-degenerate o uter hydrogen lay ers with mass ^ 
0.01 Mq at an age of ~ IGyr jPanei et al.ll2007i) . which is usually 
considered as the typical formation time of AM CVn stars. Transfer 
of this matter onto the companion upon RLOF may cause explosive 
phenomena similar to Classical N ovae, like in ordinary cataclysmic 
variables. Recently Isheni ( 1201 5h suggested that the ejected matter 
forms a common envelope in which the components may merge. 
However, to drive a firm conclusion, it is necessary to consider the 
interaction of the ejecta with the two components, having in mind, 
of course, the compactness of the system. This problem, perhaps, 
demands 3-D hydrodynamical simulations. 

Strong flashes occurring in the AM CVn systems consid¬ 
ered in the present study almost definitely can not be iden¬ 
tified with unique Helium Nova V445 Pup dAshok & Banerie3 
I 2 OO 3 L lAsholJ l2005h . since pre- outburst luminosity of the latter 
log(L/L 0 ) = 4.34 ± 0.36 dWoudt et alj 1200^ is too high 
for pre-flash CO WD in IDD. Rather, its progenitor may be a 
massive IMq) WD accreting at a rate ~ 10 ~^ Mq yr~^ 
from a massive (also ~ 1 Mq) He-star com panion dWoudt et alj 
l2009l : [^ersanti. Tornambe & Yungelsonll2014h . A possible progen- 
itor may be similar to the unique sdO-l-WD system HD 49798 
dMereghetti et alfcoogh with both unusually massive components. 
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